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Abstract. - We study the dynamical Casimir effect in the presence of a finite coherence time, 
' which is associated with a finite quality factor of the optical cavity. We use the time refraction 

model, where a fixed cavity with a modulated optical medium, replaces the empty cavity with a 
■j-J ■ vibrating mirror. Temporal coherence is described with the help of cavity quasi-mode operators. 

Asymptotic expressions for the number of photon pairs generated from vacuum are derived. 
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■ Introduction. — The dynamical Casimir effects is an instability of the electrodynam- 

ical vacuum, which has received considerable attention in recent years (see for a recent 
review pQ). For an empty optical cavity with a vibrating mirror, with no losses, it leads to 
the exponential growth of trapped photons, if the oscillating frequency of the mirror is equal 
to twice the photon frequency. Its main interest is related with the fact that it can occur 
for arbitrarily low photon energies, thus stimulating the hope of observation of quantum 
vacuum properties in small scale experiments. 

Another approach to this problem is offered by the concept of time refraction [2 , which 
is associated to the temporal change of the optical properties of a dielectric medium. The 
similarities between time refraction and the dynamical Casimir effect were recently explored 
by us [3], where the empty cavity with a vibrating mirror was replaced by a static cavity 
with a time dependent optical medium. The equivalence between the two configurations 
was demonstrated, and again, in the absence of losses, an exponential growth of the photon 
number inside the cavity was obtained. We have also discussed the influence of losses, which 
are always present in any optical cavity. These losses could be due to dissipation in the 
optical medium, or to the imperfections of the cavity mirrors. The existence of losses can 
therefore be associated with a finite quality factor Q, which on the other hand, characterizes 
the finite lifetime of the trapped photons. This will lead to a finite coherence time for the 
quantum field states, as shown here. 

Our phenomenological model for dissipation [3] had predicted that the exponential 
growth of the photon number would only take place for a photon emission rate larger than 
the photon losses. This agrees qualitatively with other models considering the dynamical 
Casimir effects in the presence of losses pQHHSj- In such models, no asymptotic value for the 
emitted number of photons would be attained for very long times. 

In contrast, recent work 0|H], has shown that an asymptotic equilibrium value could 
be derived, for low excitation rates. Such a value would result from a balance between the 
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emission process and the photon losses. However, with this alternative theoretical model, 
no time dependent result can be derived, giving us no information on the dynamics of the 
photon emission process. Even if the results were in a fair agreement with ours for a single 
cavity mode, the question remains that an asymptotic state was shown to exist in contrast 
with our previous derivation. 

Here we return to the problem of the cavity losses, and show that, by using our time 
refraction approach it is possible to derive a temporal photon emission rate which tends, for 
very long times to a constant asymptotic value. Under appropriate assumptions, it is possible 
to show that this asymptotic value coincides with the one derived by [7] , thus establishing a 
bridge between the two theoretical models and showing that the two approaches can indeed 
be reconciled into a single unified picture. The key element in this unified picture is provided 
by the existence of a finite coherence time of the photon quantum states associated to the 
cavity losses. 

In this work, we extend our time refraction model to consider the presence of losses, 
by using the concept of cavity quasi-mode, which was first developed in the frame of laser 
cavities [5], and then extended and applied by many authors [TOl-12 . This concept, and the 
corresponding field operators, are described in Section 2. Time refraction of a cavity quasi- 
mode is then discussed in Section 3, where an expression for the time dependent photon 
number creation is established and its main physical consequences are discussed. Finally, in 
Section 4, we state our conclusions. 

Cavity quasi-modes. — We start by considering a static dielectric medium in a one- 
dimensional optical cavity, which is defined by two parallel mirrors, placed at positions z = 
and z = L. The transverse mode structure is somewhat marginal to our present discussion 
and is ignored here. It can however be easily accommodated in the formalism, and it will 
be briefly discussed at the conclusions. The classical electric field describing the modes in 
the cavity can be written as 



E(z, t) = E m sm(k m z) exp(—iLu m t) + c.c. (1) 



with the mode wavenumbers and frequencies 

k m = 27r^ , u m = / fc? " C (2) 

where m is an integer, c the velocity of light in vacuum, and e(ui) is the dielectric function 
of the medium, assumed gcncrically to depend on the field frequency u. This is replaced, in 
a quantum description, by the electric field operator 

E{z,t) =Y,[E m {t) exp(ik m z) + /i.e.] (3) 

m 

with 

E m (t) = iC m \a m (t) ~ al m (t)] (4) 



where C m = \J hu m /2e(Lu m ) is an appropriate normalizing factor. It should be noticed 
that the standing wave field described by equation ([!]) results from the superposition of two 
counter-propagating photon states, as explicitly represented in equations (J2])-®- For each 
cavity mode m, the destruction and creation operators can be defined as 

a m (t) = A m exp(~-iu m t) , a}_ m = A^_ m cxp(iui m t) (5) 

Let us now assume that these modes correspond to a lossy cavity, or in other words, a cavity 
with a finite quality factor Q m . This is equivalent to introduce a finite spectral width for 
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the cavity mode, 7 ? „ = uj m /Q m . In this case, we can still use a similar mode description, 
but where the cavity modes are not really electric field modes, and can be considered as 
cavity quasi- modes. The operators a m (t) and a_ m (t) now become effective mode operators, 
also called quasi-mode operators, resulting from the superposition of an infinite number of 
field modes with nearly equal frequencies u> ~ uj m 12 . These effective operators, can then 
be stated as 

A t M _ f„. U^t <<- ^ 



Om(t) = J g m (Oa n (t,t)^ , al m (t) = J 9* m (0al m (^t)^ (6) 

where £ = (uj — uj m ) is an internal frequency variable, and <? m (0 is the spectral shape 
function to be defined later. It is obvious from here that the effective mode operators a m (t) 
and a m (t) can be considered as external or macroscopic operators, whereas the field modes 
5m(£ji) a n d a_ m (£,i) are the internal or microscopic ones. These internal operators can 
similarly be decomposed in amplitude and phase as 

a m (t) = A m (0 exp[-#(i)] , al m = Al m exp[ty(f)] (7) 

where the phase varies in time according to <f>(t) = U)t = (tJ m + £)t. Comparing this with 
equation (|5|), we conclude that 

A m = J 9 m (0A m (0eM-m§; > Al m = J g * m (Z)Al m (OeMm§; (8) 
For each internal field mode £ we can define a photon number operator, in the usual way, as 

N(0 = dl(C,t)a m (U) = 4,(0^(0 (9) 

Using equations ([7|l. we can then define a photon number operator for the effective cavity 
mode to, as 

N m (t) =/§/ ^5m(0flm(O4,(0An(O (10) 

At this point it should be noticed that the spectral shape function g m (£,) can be chosen in 
such a way as to satisfy the normalization condition 



/ i 9 .„«)i 2 ! = 1 



("I 



This particular choice will enable the effective mode operators to obey the commutation 
relation [a m , a^,] = 8 m ^ m > , given the usual commutation relations valid for the internal field 
mode operators 

MO, aL(0] = -O ( 12 ) 

as demonstrated by previous work on lossy cavities |12j . Until now, we have assumed the 
shape function <? m (£) as an arbitrary complex function, obeying the normalization condition 
(|10j) . But a natural choice, which describes a Lorentzian mode profile with half- width j m is 
given by 

(£ +7m) 1 

It is obvious that the particular case |<? m (£)| 2 — 27r<5(^) would lead to A m = A m and 
A_ m — A^_ m , and corresponds to the ideal case of a perfect cavity, valid in the limit 

7m -> 0. 

Before concluding this introductory section on the cavity quasi-modes let us consider the 
field Hamiltonian H. Following the usual definition, we can obviously write H = ^2 m H m , 
with 

H m = h I (uj m + f ) 



dti (14) 
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This allows us to derive the Heisenberg equations for the field operators, as 

^-a m (£,t) = i [a m (£,t),H m } = -i(u> m + g)a m (£,t) (15) 
at in 

with a similar expression for al m (£, t). This could also be recovered directly from equations 
©• 

Time refraction of a quasi-mode. — Let us now consider the temporal variation of 
the optical properties of the medium inside the cavity. This can be described by assuming 
a varying dielectric function e(u) m ,t). Given the fact that k m is fixed, as imposed by the 
fixed boundaries, the frequency uj m has to change in time, in order to satisfy the dispersion 
relation ([2]) for all times. This situation corresponds to a time refraction inside the cavity, 
as previously discussed by us [!1[3]. The cavity quasi- modes can now be defined as 

a m {t) = A m (t) exp[-i(f> m {t)] , al m = A ] _ m (t) exp[i0 m (i)] (16) 

where, in contrast with equation ([3]), the amplitudes are now time dependent operators, and 
the phase is determined by 

4>m(t) = J U) m (t')dt = k m C J -yjjr (17) 

where n{t) = \J e(u rn , t) is the refractive index of the non-stationary medium. As for the 
internal field modes, they can still be described by equations ([7]), but with time varying 
amplitudes A m (^,t) and JiL m (^,i), and with a new phase defined by </>^(t) = <frm{t) + £,t. 

We can now use the well known properties of time refraction [2)ll3j . As previously shown, 
in a time varying medium, the field mode operators are coupled by the evolution equations 

-TL a m = -iua m + f{t)a}_ m , ^-a-m = iuj a f _ m + f(t)a m (18) 

where lo — w TO +£, and f(t) is a real function defining the temporal evolution of the refractive 
index 

/(*) - ~ (19) 
2n at 

We can also rewrite equations (|T5)) in terms of the field amplitude operators, as 

jA m = Vi {t)Al m , jAl m = vl{t)A m (20) 

where for simplicity, the dependence of the operators on £ was not explicitly stated, and 
where we have defined 

^(t) = /(i)exp[2i^(i)] (21) 
Integration of equations (|2"0")) leads to the following solution 

A m (t) = cosh[r(e,t)]i m (0) + sinh[r(e,t)]i t _ r „(0)e^ 
Al m (t) = cosh[r(C,t)]iL m (0) + sinh[r(e,t)]i m (0)e-^ (22) 

where we have introduced the squeezing function 

r(£,t)=[ \vS')\dt' = J f(t')e 2lit 'e^^^-^dt' (23) 

The phase 5^ is defined by v^{t) = \v^{t) \ exp(i(5^), and will not contribute to the mean photon 
number, to be considered below. We can therefore take it as 5^ = 0. An alternative and 
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<N m (t) >=2 7m / ^%V J £ (25) 



equivalent option, would be to incorporate this phase in the definition of the operator A^_ m . 
We can now calculate, for each mode U) or £, the mean photon number, as < 7V m (£, t) >=< 
0|aTj a TO |0 >, where |0 > is the unperturbed vacuum state. This leads to 

<N m (U) >=sinh 2 [r(£,t)] (24) 

It should be noticed that the same number of photons will be emitted in the opposite 
direction, < N^ m (£,t) >=< N m (£,t) >, in order to guarantee that the total momentum 
of the electromagnetic vacuum is conserved and remains equal to zero. This expression 
therefore states the number of photon pair emitted from vacuum. According to equation 
(ITU)) , this then leads to the following result for the mean photon number associated with the 
cavity quasi-mode m 

sinh 2 [r(£,i)] d£ 

where a Loremtzian profile was assumed. This is our man result, which has important 
physical consequences to be discussed next. In particular, it will lead to the occurrence of 
a saturation of the dynamical Casimir instability, due to the existence of a finite coherence 
time for the cavity fields, charaterized by lj^j m . 

Let us first assume a sinusoidal perturbation of the refractive index of the medium, such 
that 

n(t) — no + Sn(t) , Sn(t) — noesm(0t) (26) 

where O is the driving frequency, and the amplitude perturbation is assumed very small, 
eCl. As discussed previously by us [3] , this is optically equivalent to an empty cavity with 
an oscillating length, L(t) = L + SL(t), where SL(t) = L 5n(t)/n - In this case, we can 
write 

/(t)~^cos(fit) (27) 

ft 

4>m{t) = I u m (t')dt' ~ uj t + cos(fii) (28) 
where we have used ujq = k m c/n Q . This leads to a phase factor of the form 

oo 

e 2^ m (t) = i l J l {a)e^ +l ^ t (29) 

/ — — oo 

where a = 2eojo/n, and J/ are Bessel functions of order I. Replacing this in the expression 
for the squeezing function, we get 

r(U) = ^ y)i'Ji(a) f W w+t ' +e iu ~ t '] e 2 ^' dt' (30) 
4 i J 1 J 

where 

cj± = 2w + (l± 1)0, (31) 

In order to maximize the value of r(£,i), we have to choose the driving frequency O such 
that one of the frequencies u>± is made equal to zero. This corresponds to 

= O t = (32) 

which implies that 

a = ai = -e{l±l) (33) 

We generally expect to have a«l, even for large values of |Z| 3> 1. This allows us to use the 
approximate expression for the Bessel functions, J/(a) ~ a l /2 l U. It becomes apparent that 
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the dominant term in equation ((301 will correspond to I — 0, although this is not convenient 
experimentally because it implies a large driving frequency fl = 2u> . The squeezing function 
then becomes 



Jo 



(34) 



with 2/0 = Let us first consider the case |r(£,t)| <C 1. This allows us to simplify 

equation (]25[) . which can now be apprimately written as 



< N m (t) >= v$ 



2 i d£ 2 7 „ 



2^ {e+ii) 



t n 2 



This can also be stated as 



Noting that 



<JV m (t) >=vl fdt 1 f dt" f ^ 2% 
Jo Jo 



dio 



27T (e+l?n) 



™ e 2if(t'+t") 



(w2 + 72) 27r =^ eX P(-7l*D 



integration over £ can easily be performed, leading to 

< N m (t) >= vl [ dt' [ dt" f e -^l*'+*"l 
Jo Jo J 



e-i mt 'dt' 



We can finally write 



<N m (t) > = ^-(l- e -^*) : 



(35) 



(36) 



(37) 



(38) 



(39) 



We should notice here that, for f>T, where r = l/7 m can be define as the field coherence 
time for the cavity mode m, the number of photon pairs determined by this expression 
attains the asymptotic value of 



< iV m (00) >= vl T 2 



(40) 



This eactly coincides with the steady state result derived from a completely different method 
by Lambrecht et al [7], for a single cavity mode, thus showing that our time refraction ap- 
proach to the dynamical Casimir effect is indeed capable to describe the temporal evolution 
of the number of photon pairs created from vacuum, and to establish an asymptotic state 
for a lossy cavity. The possible existence of different modes contributing to the effect will 
be discussed below. On the other end, equation (|3"9l also shows that, in the oposite limit 
of short times, t < r, we get < N m (t) >= t 2 , which agrees with our previous result in 
the same limit [3] . We can therefore state that the existence of an asymptotic steady state 
is provided by a finite decoherence time r for the field in the cavity. 

The above discussion is only valid for vqt -C 1, compatible with the assumption of 
|r(£,t)| <C 1. which led to equation (|3"9"|) . But, in order to observe photon creation, and 
to get < N m (oo) > > 1, we should be able to consider a more general case. Let us then 
determine the number of photons < N m (t) > for an arbitrary squeezing factor. 

In the general case, where |r(£, t)\ can take values much larger than one, we can develop 
the sinh 2 function in a power series, and for each term of this series, perform the integration 
over £ in the way described above. After taking the sum of the resulting series, the final 
result is 

r vo 



< N m (t) >= sinlr 



(1 



(41) 



We can see here that the initial exponential growth associated with the dynamical instability 
is limited by a coherence time. This result is therefore qualitatively different from previous 
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Fig. 1: Number of photon pairs created from vacuum: (a) as determined by the phenomenological 
model < N m (t) >= smh 2 (uot) exp(— j m t), where exponential growth always occurs for uo > j m ; 
(b) as determined by equation <41D . which shows the existence of a asymptotic saturation level. 

estimates based on phenomenological models for the cavity losses, which always predicted 
exponential growth, as long as vo > 7m- This is illustrated in the figure. Here again the 
same condition is required to excite a large number of photons, but this number tends to 
saturate for long times, f>r. The asymptotic value is now given by 



which generalizes the asympotic expression derived by [7]. 

Conclusions. — In this work we have considered the emission of photon pairs due to 
the occurrence of time refraction in a lossy cavity with an oscillating dielectric medium. 
This configuration is equivalent to that of a dynamical Casimir effect in a lossy cavity with 
an oscillating length. This extends our previous work [3] to the case where the coherence 
time of the electromagnetic field in the cavity is finite, as imposed by a finite lifetime of the 
confined photons. Such an extension was based on the use of cavity quasi-mode or effective 
mode operators. 

We have shown that our time refraction approach to the dynamical Casimir effect is 
able to describe the temporal evolution of the photon number emission from vacuum. In 
particular, it was shown that the existence of a finite coherence time leads to a photon 
number saturation. Our result for this saturation value is valid for arbitrarily large squeezing 
parameters, and reduces to the value previously found from a completely different method [7] 
for a single cavity mode, for the particular case of weak squeezing. However, in contrast with 
[7] , we only consider processes associated with photon pair production with equal frequencies 
u = £1/2, and we exclude processes associated with unbalanced photon frequencies lu + lu' = 
il. The reason is that these unbalanced emission processes violate the conservation of 
total vacuum momentum. Of course, we can always excite more than one cavity mode. 
First of all, we have modes with equal axial wavenumber, k m , but different transverse 
structure. In general terms, this transverse optical structure can be decomposed into a sum 




(42) 
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of orthogonal Guass-Hermite or Guass-Laguerre states [14]. These modes have a slightly 
different dispersion relation, and a bunch of nearest cavity modes can eventually be excited 
inside a short frequency band, as imposed by the oscillation of the optical medium. But the 
losses for the higher transverse modes tend also to be larger. On the other hand, different 
axial modes k m can also be excited, as defined by the condition (|32[) . But, under plausible 
experimental conditions, only one cavity mode will dominate, as considered here. 
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